MASSACHUSETTS INSTITUTE OF TECHNOLOGY 
PROJECT MAC 



Artificial Intelligence MAC-H-371 

Memo. So. 159 April 196( 






NUMERICAL SOLUTION 

or 

ELLIPTIC BOUNDAKV VALUE PROBLEMS 

BY 

SPLINE FUNCTIONS 



Jaytint M« Shnh 



A numerical method for solving linear, two~d1mo.n&jonat elliptic boundary 
value problems is presented* The method is essentially the Rit?. prmrodurt 

tAich use* polynomial spline functions to approximate the exact solution. 
The spline functions are constructed by defining a polynomial function 
over each of a set of disjoint subdomains.and imposing certain compatibility 
conditions along common boundaries between subdomains, The main advantage 
of trie method is that it does not even require the continuity of the spline 
functions across chr boundaries between subdoraalns. Therefore it is easy 
Co construct classes of spline functions which will produce any specified 
rate of convergence* 






i> INTRODUCTION 

During the last few years polynomial splines have been ust*d to obtain 

numerical solutions of specific elliptic boundary value problems [2,3,41 

This paper presents a general formulation of this type for solving tvo- 

dimensional elliptic boundary value problems* The method is essentially 

the Ritfc procedure applied to a finite ditnonsion.il apace of polynomial 

spline functions. The formulation can be extended to higher dimensional 
cases and to simultaneous systems. 

I 

Ue denote by C a bounded open domain in the Euclidean plane with boundary 
*G and closure C. For a function ugC (G) we denote by I | u I I U r ^ c maximum 



1/2 



norm max max ID u(x) I and denote by I lull, - the norm 



J 2-f |D a u(x)| 2 dx . 



a 



Here a = (a., tO is the multiple index of length |a| * c^ + c* 2 and D is the 

corresponding partial derivative. Wo denote by |u | , _ the semi-norm 

2-1 ID u| 2 dx) . U't L = X, b WD 3 be a real strongly elliptic 
|a]-i ' / la( s 2m a 

operator of order 2m defined In G. Suppose that b„gC (C) and G is of class 
C , Let the bilinear form associated with the operator L be given by 



|a|*m -n 
IP 



2 * 






- £ (-1)1*1 D a (bJ). 






where u fV »jC (G>. Then L - £^ ('O^ I D (b .IT). Wo assume tliat B(**0 
satisfies the Garding's inequality in the following form: 



(1.2) | B <v»q£|* c M m , G for a11 q> e C*(G') 



where G r Is any open subset of G. 

For a given f € C (G) , let ucCq(G) be the infinitely differentiate solution 
satisfying the equation Lu ■ E end satisfying zero Dirichlet data on the 
boundary. Then 



(1.3) 



" 



, u) = f ^dx ^CjCG) 



2. TWO- DIMENSIONAL SPLINE FUNCTIONS 

Subdivide the domain G by a grid h into a finite number of subdoaains 

* ■ 

such that each subdomain is a polygon vith infinitely dif ferentiable curvilinear 
qides* Wc call any such subdomain a cell. Let (S f ) be the set of grid lines 
forming the sides o£ colts. Wc require that the interior of any grid line S 

■ 



must not contain a vertex, L*t sup (Length of 5.) " h. Wo denote by [G.] the 

disjoint -rollectlon of Interior* of the colls and let G fl = \j G . 

i 

Assume that ejeh grid line is given by an equation of the form 

■ 

(1.4) x - f ± ( ?l ) 

y - 8i Uj) 

where £ is the distance along the grid line S ractirtured from on*? of the 
vertices which It connects. The unit tangent vector c. - (t- *t ) at a point on 
the grid line and the unit normal n, ■ (-t ? ,t ) are defined by the equation: 

df./de. 
d.5) t = r ^— * j-^ 

r df i % 



dg,/^ 



r i 3T> + ( dT' 



If u 4 «C X (G. )> let D*J* u (l ' denote the r-th derivative of u* 1 * in the 
i n,r 



■ 

4 



direction of the normal n<j defined tf S. is one of the sides of G . For every 
cell, w define the following quantities: 



<>•<> ■£ • jk I 4 « •« ««, 



defined for all k anc r 
whenever 

Ostsm-l 
and (kk^p+r-(ra-l) 



where p id samo fixed integer ^0 and Independent of |* , Also if w is a vertex 
of a celt G t , dofioo a* ,w - D a u tlJ (w) for all a such that 0^|a| s ra-2, 

Wc define a ciass^(i},u) of functions with domain G as follows! Let 
lK/»yj) denote the local coordinates of a cell G, , obtained by translation 
of the origin to a point in G., If V \>- ■ tnen u = u restricted to 0. 

is given by an equation of the fern 



(1.7) u (i) - a oo + a u x L + . tfl y. + ... + « y" 1 



such that the following compatibility conditions arc satisfied: 
(i) For all fc and * and whenever defined 



(l - 8) \r " V j if S j -a^Wv 



and uJ-J. = If S - aGfaGj 



(ii) At a vertex u, for all laU ra-2 



(1.9) u*' w - uj ,w if w belongi Co both G t and C 1 , 



i v 

and u ' "0 if the vg^G 



Ko also require that 
(1. 10) n a n+ p 

Mow we impose the following conditions on the class of admissible grids: 
First we Lntroduceron-dimensional coordinates X. • x./h, Y, ■ y./h, *. ■ £. . /h 
and quantities A, k - a h Jtk ,D ££ - u^ andU*' w - u* ,u hl a l so that 



(I.ID u<« =A Qa+ A lo X 1 + ...A y" 1 

1 
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d.iaj $ - | = J »^S. %' w - frHtfA 



Far a given cell G. f let V bu the vector space spanned by the quantities 
U.j! and IT* and let A, be the vector space spanned by the coefficients A,,, 
Then we have a linear transformation 






(1.13) T £ : A t - V t 



\ 



We choose n sufficiently l«rfl? so thaC T, is onto, Hpw let K. be the kernel 
of T. and let A, be the conplcment of K . * Then T induces ,m isomorphi 



n 



(l.M) T t : A^Vj 



We require that ; (a) Values of n. are bounded uniformly with respect to all 
admissible grids and the absolute values of the eigenvalues of T. are bounded 
away from zero uniformly with respect to all admissible grids, (b) There exiats 
a positive constant : such that for any admissible grid, 
lenRth of S, 



inf — — r L a 5 

t h 



I 



(c) There exists a positive constant & such that for every admissible grid, 

and for e*/ery ceil, there exist* a point (x <y ) in that cell such that the 

square {(x,y)||x-x l-rj ? h, Jy-y [sfi-hj is contained in the ceil, (d) Modulo 

translation* rotation and scale factor, grid lines S, arc restricted to a 
finite number of shaped, (e) There exists a constant C such that the number 
of cells that any straight line through G intersect* is leas than ^q/' 1 - 

If the cells are triangular or rectangular in shape with straight sides 



and if they satisfy conditions (b) and (c), then it is possible to choose a 
single integer n for n. in Eq. (1*7) such that condition (a) is satisfied for 
all such cells. To show this for triangular cells, consider a cell G in the 
form of an isosceles right-angled triangle with hypothenuae h units lonfl and 
with origin of the local coordinate system at one of the vertices. Since the 
transformation T. for this particular cell doss not depend upon h, we can find 
an integer n for n. in Eq. (1.7) such that condition (a) is satisfied for all 
values of h. flow consider any other triangular cell G f which satisfies condi- 
tions (b) and (c). Set n , ■ n. Without loss of generality we can assume that 
G, has the origin of the local coordinate system at one of its vertices. There 



exists a notrstngular linear transformation S mapping C- ; onto G, , * Nov 
specification of the quantities u * v and u. ^ {i£ defined) Implies #pecl flec- 



tion of all quantities of the form 



k+1 



^ *J n,r t,r" i 



O^r+r'^m-1 

r'>0 

Oiik^+tr+r^-fm-l) 



where D^ , u^ ' denotes the r*-th derivative of u 4 J in the direction of the 
t it 

tangent, defined if S. is one of the aides of C . * 
Therefore S Induces en isomorphism 



a. i7> 



S 2 :V i~V 



Also 5 induces an isomorphism 



(1.18) 



S. : A - A £l 



Therefore we have a commutative diagram 



(1.19) 







V i > ° 



1 ;{ ,1 



-> V > "i'-=—> V > 



1 T. ' 

1 



u 






This , In Cum, implies a cocmutfitivo diagram 



(1-20) A t ^ V ± 



I *■" 



i' 



■ 






Therefore T, , is nonsingular. Now we claim that condition (a) is satisfied- 
For, if not, there exists a sequence of triangular cells [G. j such that the 
corresponding sequence of the lowest absolute eigenvalues of the transforma- 
tions T tends to zero. But any triangular cell is specified by six quantities, 

K 

namely the non-dimensional coordinates of its vertices. Therefore the sequence 
of triangles can be embedded as a closed subspace lying in the closed annulue 
In JR , lying between the unit sphere and the sphere of radius g-. Since the 
annulus is compact, the sequence G. has an accumulation point. By continuity, 

the triangular cell (V corresponding to the accumulation point satisfies con- 

* 
dicions (b) and (c) and has singular T fco which contradicts (i.20). 
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In the sane fashion, we can choose a single Integer in Eq. (1.7) whixh 
will satisfy condition (a) for all rectangular colls provided that condit.ons 
(b) and (c) are satisfied. 

For practical applications, it suffices to segment a draoin using tri- 
angular and rectangular cells aa well as a finite number of special shapes to 



handle a curved boundary. 



The vector space *Jf(h,p) is constructed as follows. For each coll G, » we 
have an isomorphism A ^ V ©K . Then a function ■ ^. r b, is defined by 
specifying the quantities 



kr Tcr fc t 



4-4r -■£ '"-••• 



* 



for each grid line not 
in gG 



for each vertex not 
in 3 G 



and an elepent of K for each cell G 



Ue require that 



u^ r = tf S is in aG 

i v 

and u ■ » if w is in gG 
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■ 






We extend the (semi) norms J- 1 » 1 1- 1 1 > I 1 1- 1 I ami l ' ,c bilinear EVrw B(*r) to 

pit'*?} by restricting the integrals to C n . 

y 

Let V be a vector space spanned by the quantities u~ > u (S. and w not in 



1 



5G) with the basis {e^. Let K= ®K t with basis [k ] . Lot W V g K. Ihon 



we have an isomorphism I t fO >p) ^ W r The images of e f and k, under I foro 
a basis for ^( ,p). Denote this basis by £u ] 



3, FOBMLTLATION 

■ 

We seek a solution u e ^(i ,p) such thnt the equation 



(<p>"> ■ J r cp f 



(1.21) BU,u) = I ,'fdx 



is satisfied for all &&;*{ *p) • W- show that Eq. (1.21) has a unique solution 
d that D Q (u-u) - 0(h p+ T' for Os|a|*m-l . 



an 



For a numerical solution of Eq. (1.21) we reduce it to an equivalent matrix 
equation as follows* We want to find v £ U such that u ■ I~ (v) satisfies 
Eq. (1.21). Carrying out Che Integrations, B(m»u) ■ (v,B v) where v ■ Kg), 



12 



( t ) is the loner product and B is a matrix (&,•) where B « B(u ,u )• Also 
/cfdx ■ (v»FJ where F € V and Its i-th component is given by /u.fdx. Therefore* 

Eq. (1*2L) is equivalent to the matrix equation: 

■ 

- 

(1.22) B v - F 



■ 



If within each cell, vc approximate the coefficients b*- and the function 
f by polynomials use of quadrature formulas for integration of polynomials over 
eimplcxo5[5] simplifies the computation of B- The computation is carried out 

by computing B(u >u ) cell by cell. 

■ ■ 

4 - COWVERGEECE 

In the following C 1 , C, T t . .. denote positive numbers , the choice of which 

does not depend uponh - 

Lemma 1: Let u e Y(h,p). Let S, be a grid line contained in jG,ngG f ,snd 

let u ■ u restricted to G and u * u restricted to G f , . Then for lal^m-i 



(1.23) 



«*£*, ftA •?> V<x."h" n " 1 " |u| < hU^I"!.^ 






u 



<1.»> 



If S. - ^G-reC, then for |a|^m-l 



oJ&, n/w>W^M l»U ( 



Proof: Suppose that S, connects vertices w_ and u^. Without loss of geneva 11 ty 

we con assume that w. Is the origin and that the x-axis passes through w^. Let 

X, be the length of S 4 * Now consider a subspace d °f the vector space A. 

spanned by (A QO , A u A fll * --- * 0taH>I 3 and a subspaceb> of V t spanned by 

flf 4,w ll, , ft . (U 1 ^!. , - and U 1 -* , . (If ra - L, then U is one- 
t a J a 4**2 ' L a JJajsm-Z o,m-l * o 

dimensional, spanned by U Q >- Let J : Q. - Aj be the inclusion map and 

p ; V . -. (j» be the projection map. Then It Is easy to check that the map 

?, : Ci -foi given by the composition pot oj is an isomorphism and therefore 

°l cl . Moreover since grid lines are restricted to a finite number of 

o i 

♦. 
shapes and \./h € [j).*!]* It fallows that the absolute values of th? eigen- 

values of T . are bounded away froa zero uniformly with repect to n . 

Let I*. be the complement of (> o in V and Q ( be the complement of Q o 
In A. - Lecd/ »<%■# 6* f l>i be the corresponding spaces defined for the cell G £r . 
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Ut 



u s T*(k,p) and lot w c W and l; c K bo fcho wctura corresponding to u- 
Let the projection of w onto V, be v and let v » v © v 1 c byD t^ i*e« v^ | >o 
and v.cti- LeC a - T " l (v ) and let the unique polynomial corresponding to a^ 



be F = p + p„ x. + .•■ + P ^Y** . Similarly, the corresponding vector 



v^ € t^ determines a unique polynomial P' = p QQ + pJ Q x i **^o l m-l Y i * BUt 



v = v 1 - Therefore P = P' . 
o o 



LcC =" 1 



d = Tj'CCodWj)) and d - (d^,) e C^ffiC^ «*erc d « rt and d^-l^ . 
Let T * » T restricted to a , q - *>,l. Then poT Q + poTj is the zero map and 
poT* - Ti * Therefore d Q - - T J l (tJ (dj)) »*ich inplieB that if d Q - (V QO , 



no v*-i } and d i = c vc •••• W' *« 1^1*2 ^ "V £or 



|ot|<». 

We define the following norms: 



Por q - < %Q , q l0 , --.. 1,^*1^ «■«« che no ™ I M I^J Kj I 1 

Box r - (r 00 , r 10 , ....r^.Oi defioe the norm | |r||- ■« ^ || 

<*«,o "• W**' de£tne the M ™ IIMl-J^hjl 



ir s 



■ 



■ 



■ 

* 





v 






r 


' 




* 

■ 
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. 












■ 




.. 


■ 










, 
























, 











For r^fl ■ If <MK| define Che norm | | z \ | - \\r\\ + ||i[|- 

Then l^ilrfjIldjU. Similarly | |d;| [rftj 1 |dj || for the cell flj, - Now 

■ 

T.[((a + d ) © d.)) ■ v. Lot the projection of k onto K ± be the vector 






■ ■ 



v - (k i . >. k ) and let d. e K " £, * nd d ! © k ! " Qj- Thon 
1 njiO 0|Tii l i j. i** 






•i 



■ 

I 

■ I 



■ " C "- ^^•■•^.m-rt.m-!"? ^.X.X*-^,^;..!^ 



■ 









! 



u U). u (» , )= (Yoo . ¥;o H...^ Yo ^ 1 . v ; i0 _ l )^- l ^k^ o ^ rai0 .k; i0 -v; i 






Therefor* sup lD r ,u (I) (x,y)-D nU ti W*r-jg7 ( I 1 1 £ | | + H^, ||) 

(x,y) e Sj a hi I 






I 

Sow we claim that 

- 



- 






(1.25) llMl SC 4 h,i ' l 1 U U i 



■ 



To prove this observe that 



■ 



. 



1 H.. 0l -[/^(^-) *i«J 



















- 


















■ 








■ 














■ 
















1 



■ 



■ 

J 

' - ■ ' ' - ' ■■■ , - ■ ■ ' 



■ 

16 

r v* ( A*& 



,£ 



16 

■ 

- - ■ ■ 

2 



The quadratic for* [ , 1=2-* ( & a , * 2 ) ^l^l 



is positive definite on the linear space CVgKa - Therefore 

inf U2 

OHttXj&t -ferj " \ I* - \ U a ^nctloa of 1 and \) - Suppose lnf X ( - 0. 

Then there exists a sequence or cells {<y such that the corresponding 

sequence \ - 0. Now each G k contains a point (x^y^ and the square 

■ - <* ■ {(x^y>ll*-\l* n v! y " y k^ ::h6 2 , * ut ' 1 ' v denot0 the MUncar £om [ ■ ' 

restricted to the square o^. Then there exists a sequence of points < V V * 

- 

l^lsl-fij.plJail-Sj and a corre9 P° ndin S "1 uen « i^} » uch chflC ll'ltH" 1 and 

■ 

k^l'i 0, By compactness, there exists a subsequence {\ >\ *\ ) c 
verging to (^ ^ ^ > in iR 2 ©^^ «•* that ||^|| - I and tZ^J' ■ 
which is .1 contradiction of the posltlve-definiteness of [ , 1'- This proves 



' 



on- 









the claisw Inequality (U23) follows. The proof of Inequality (1.24) is similar. 
Lemma 2: Let Uc3fchfP}> Then 



L 



<'•") lll»IIL-i,c o , | | u |„, 0o 



■ 



, 



. 









. 






. 17 

■ 

Pt „„„ WUhou. loaa of ..nora.ity « m, ..«- *.. none of m *U •'-< 
ls p,ralUl to the coordinate a,... Mi- - ft-** "' ""^ * U "" n8 



. 



„, . „ „„ C, and »■ - «lae>*er e . Tor any Urn , - ?„ int«acctio„ Q . f to* 
point, (»..,„> and <»„.,„> in S « •* « M *«*» k"« ""^"' t " 
<,.,„> and „in B in be b.t«« <« a .y > - W" *« ■ «-' « " ( °i 'V ■ 

| 

lal -10-1. define a function 



- 



otherwise 



- 



Then J a Ct) U fUc«l« continuous ,n<- differentiae over open se^nt 



. 



I 



Define J (t t +) ° ^i» J Q (« 



t-t t >o 



and J a (t t -> -Li«J a <t> 
t-t ± 
C-t t <0 






Let ^ - ^ for C in one of the open sclents listed above. Then 



. 



• 






* r. .■*■*'. 

■ 
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J (t -> - J tf (t 1+ ) + J jy T )4 T fi«i-0 i l,....k-l. .". For any point 






tett.-t^,), 









• . 






1 

L 






I . 






' 



- 






jj ! ■ ' I 

■ 



t 



J- 

- 






■ 



■ . ' 



*2 






! 

i 









* 2 






■ 



by Cauehy-Schwarz Inequality- From Inequalities (1.23) and (1.24), 

2 



■ 



■ 



' 



From Inequality (1.25), l J a <T, l * "i l U 'm,G If 

h It 



>,£*,.,* 



. 



h* "'"k 









... . . 

■ 



- i ■ 

■ . 

, ■ 

. * 

^V*^ " xJ* yJ € G h Therefore 

P B p V K 

. - 
G 
Condition (e) oq the clas? of admissible grid? implies that i £ o and 



compactness of G Implies that t£C_ 

2 
Therefore | J (c> | 2 * C 8 Iu| 2 


- 

- 

■ 




h »,G A 


- 


«d IIIVHInG " M|J Q MI°.ie s -£ l"l 


for |a| " bi-1 


h* m > C o 


. 


flv aucceiksLve Integrations In a similar fashion^ 











IIIV'Ifo.G * -f H-.B f " I"' * - 1 ' 

The lemma follows. 



■ 



■ ■ . ♦ ■ 

Wo state without proof the following special form of the Green's 

■ 

■ - 

formula. The proof is similar to the one given in reference [l] 



■ 



: 












i - 

■ 

■ 

■ 

■ : 
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Lemma 3. There cxlat linear differential operators W R (x,y.D) tor (x.y^Ud^ 
with coefficient* bounded uniformly «±th mp<et to all <x,y) ami H ul,, .u.-h 
that H Is of order £k and for all ueC = (G) and {s^yiv^?) w 



».« | j. - B V uHE g J J,(» aiJ 



W-D nJ ^W^l-j^ 



-i r .<i k > 



iRi^'W 



s w «ac 



where the grid line S = 5G 03C! ' or 3GD 3G 

k l k *k k 



(1.26) jf »Lu- B(#,u)| < C U hP+1 t#l m ,C o for all * cy(Sp) a nd all^ 



° 



Lcnaaa *l. For a given a e C (G) 

P+li 

Proof: If o - (Qj, Oj). let al - e^! a,I and 

(x, yf = x y Then by Taylors' theorem 

* " E 7t * U <V 'a' U ~ X o )Ql<y " V' + Rtx » y) 
|d [tfp+m 



21 



where | | |lP R(x f y) | | | ^C 10 <i P+ID+l *H and 



d = distance between points (x,y) and (x , y ). 

Nov any grid line S. U given by the equations x - f^ U fe ) » y ° RfcC^)- 
Let (x , y ) be a point in S, - Since f. and g. are C functions and since 
the number of shapes of grid lines is finite, It is possible to choose 



C ?n such that 



(1.29) D ttJ u - j a^ + R n (£ R ) 



r-0 

p4TO-j+l 



where MIWIH^* C 10 £ k 



Substituting U#29) in Eq» (1.27) 



■ . 

, ■ 

■ 

■ 

■ ' 

/ s^G j-0 r-0 . k 

■ 

- 



• 



ra-1 p»j+J-n. ± k 



- - 



■ ■ 
■ 

i 



s.eaG j-tf r-o 







■ 




Jvj ^"i.jttfc*}] 




8 k 






: 



■.. 



. /e ill , r f ^J + 2-n> 



w. ie n lUi^ <«fc> ' » l - k * C 10«k 



Since | satisfies the compatibility conditions (1.8), 



■ 






1 k k 1 i k 






if s. - 3C, 0« 



- i * 

■ 

By Cauchy- Schwa rz Inequality, 












■ 






■ 
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!j;(<v/ k -v/r 



t^y-w'^]^^] 



V vv ■'k 



£ C 



. ■- 1 ****{ s,i ii <% + <° H ' a* 4 * w *"^ 



< U8 in8 Inequalities (1.23), (1.24)) 



,< 2C 



l = J o[*"- 1 - i ^l'llo ][- W2 - mJ ' 



(area of G) "} 



< c h 2(p+l> UI 2 



Slollarly 



{S/v^Sj^} 3 



V M S H 



' C ll h 



2(p+l), A ,2 



'm,G 



There 



fore J *Lu - B(*,u)| S C H h 1 * \p\ m$ , 



' X 






.. 



■ 



24 



'-. " 



■ . : 

■ 
■ 

■ 



Theorem 1: Equation (1.21) has a unique solution u e >(*%p) and 

if 



(1.31) |||*-»HU-l t G * C,.h^ 



26' 



\ :■ 

■ 

Proof: Suppose for q,,u £ ^^*p), B(p,u) = for all v * 






- ; .-. 0- |wm>i*|££, 0o 

■ 

IIMIU-i.fi * ° by Im ^ Bli ^ (I " 26> 

u = 

Eq. (1-21) hat a unique u e ^<^P>* 
Using Taylor 1 * theorem, within each cell, we can write the exact solu- 
tion at u - , u(x l ,y 1 ) + Hx^yJ where 'u^ .y^ is a polynomial of degree 





» 


and 






















1 








■ 










. 


. 




- 








» 










■ 

■ 






















1 

■ ■ ■ 






* 










■ 








■ 



„ lM d HVfri^Ha^ 



l+fl,-|Q| ,__ ,„, _ 






-. C w h™i I"! for |u|iH t 



Then |"UJJ | s C 14 h°l +l .nd p$f\ rf^ft* 






- 









' " 
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i kl ° ) all j.»,k,t,a * 



-1 



'"go' -••• "0,1^ 



ni+i 



Then | Aij | iC^ ||d|l ^ h 



ut "ucx^) ■- 00 +* 10 x i + ... s, n £* 
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